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"<C ' Abstract 

CN| ' The in- medium elastic scattering qq — > qq,qq — > qq and qq — > qq is calculated 

within the two-flavor Polyakov-loop-extended Nambu-Jona-Lasinio model. The in- 

^$3 ! tegral and differential quark-quark scattering, its energy and temperature depen- 

*~+ dence are considered and their flavor dependence is emphasized. The comparison 

with results of other approaches is presented. The consideration is implemented to 

the case of quark-pion scattering characterizing the interaction of quark and hadron 

phases in a kinetic multiphase treatment, and the first estimate of quark-pion cross 

sections is given. A possible application of the obtained results to heavy ion collisions 

J> . is shortly discussed. 
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3T ■ 1 Introduction 

CO 

To describe high-energy nuclear physics, the knowledge of the in-medium be- 
K/ • havior of quasiparticles such as quarks, gluons, mesons and baryons including 

their antiparticles is crucial. The quantum chromodynamics (QCD) seems to 
be the best tool to proceed to this study. Nevertheless, it is well known that 
the direct implementation of the QCD Lagrangian is not feasible in this state 
except for some particular cases. In order to avoid the QCD difficulties, some 
effective models were developed. 

In this respect, the low-energy particle sector is well described by effective 
chiral theories of QCD, the Nambu and Jona-Lasinio (NJL) model [1]. The ad- 
vantage of this model is that it can be studied in the entire temperature range. 
The NJL model also offers a simple intuitive view of chiral symmetry break- 
down and restoration via the realization of the quark-antiquark pairing similar 
to the BCS theory of superconductivity. However, a simple point-like interac- 
tion form of the model does not ensure its renomalizability, and a cutoff scale 
A must be introduced in the theory. The impossibility to treat confinement- 
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deconfinement phase transition and the absence of gluons are other important 
defects of the NJL model. 

To eliminate partially this defect, it has recently been proposed to couple 
the quarks to a Polyakov loop [2,3] as a mechanism that could simulate the 
confinement, even if the model does not consider the color degrees of freedom 
as done in QCD. This realized approach is called the Polyakov-Nambu-Jona- 
Lasinio model (PNJL) [4-7]. Some recent results show that this approach 
provides some advantages [4,7]. In particular, the extended model allows one 
to correctly reproduce lattice data of QCD thermodynamics [5,8] as well as 
to improve the NJL model at low temperature due to the suppression of the 
contribution of colored states. In addition, the PNJL model is more efficient 
for describing the restoration of the chiral symmetry by a rapid decrease in 
the effective masses of the quarks [8]. 

Allowability of quark-gluon degrees of freedom along with hadronic ones means 
that the model for heavy-ion collisions should be multiphase in nature and 
include possible phase transitions between these phases. Generally, this com- 
plicated situation can be described in terms of hydrodynamics or kinetics 
which have their own advantages and disadvantages. The use of kinetics for 
the quark-gluon phase needs knowledge of in-medium cross sections for its 
constituents. In A MultiPhase Transport (AMPT) model [9] this phase is de- 
scribed in terms of the parton cascade with the partonic elastic cross sections 
estimated within the perturbative QCD (pQCD). The effect of surrounding 
matter was roughly included by introducing the effective Debye mass. 

The microscopic quark dynamics is studied in a more elaborated way in the 
Parton Hadron String Dynamics (PHSD) model [10-12] where the plasma 
evolution is solved by a Kadanoff-Baym type equation. Here the potentials 
between the plasma constituents are chosen in such a way that the model 
equation of state is consistent with lattice calculations. Cross sections are 
derived from the spacelike part of the interaction and are employed for the 
scattering interactions among the plasma constituents. In this model, gluons 
as well as quarks acquire a large mass when approaching the phase transition. 
Therefore, the prehadrons which are created in the phase transition are rather 
heavy. Another model which allows for these studies is a gluonic cascade re- 
alized in the Boltzmann Approach to Multi Parton Scattering (BAMPS) [13]. 
The gluon emission and interaction during the expansion stage of the QGP 
move the system towards equilibrium. 

The quantum molecular dynamics of the expanding q/q plasma has been pro- 
posed recently [14]. Properties of quarks as well as elastic scattering cross 
sections were calculated within the three flavor NJL model. 

All the kinetic multiphase models mentioned above should describe a tran- 



sition from one to another phase: from quarks-gluons to hadrons, in our 
case. This smooth transition is simulated by a possible coalescence of quark- 
antiquark or three quarks, being close to each other in coordinate and momen- 
tum space, to form a meson or baryon thereby creating a mixed parton-hadron 
phase which is a typical feature of the crossover phase transition. However, 
none of the models takes into account a possible interaction between different 
phases. This interaction may change the order of phase transition and result 
in observable effects (for example, in dilepton production). 

In this work, we want to make a step towards account for the interaction 
between quark-gluon and hadronic phase. Basing on the PNJL model, we give 
here the first estimate for interaction of quarks/antiquarks with pions which 
is expected a dominant component in this type of interactions. 

Thus, the purpose of this paper is to calculate quark-quark, quark-antiquark 
and antiquark-antiquark cross sections and generalize this approach to the case 
of quark-pion scattering. The consideration is based on the chiral two-flavor 
PNJL model the key points of which are remained in the next Section II. In 
Section III and IV, the main equations are given for different channels quark- 
(anti) quark elastic scattering and are generalized to the quark- hadron case in 
Section IV. Their numerical results fo qq and qH processes are presented and 
discussed in Section V. We conclude the obtained results in Section VI. 



2 The PNJL model used 



The deconfinement in a pure SUf(2) gauge theory can be simulated by in- 
troducing a complex Polyakov loop field. The two-flavor PNJL model is used 
with the following Lagrangian [15-17] : 



£pnjl = q (ij^ -m )q + G [(qq) 2 + (qi^fq) 2 } - U (<&L4], $L4]; T) (1) 

where scalar and pseudoscalar interactions are taken into account, G is the 
coupling constant, f is the Pauli matrix in the flavor space, q and q are the 
quark fields (color and flavor indices are suppressed), m is the diagonal matrix 
of the current quark mass, rfiQ = diag (m°,m^) and m° = m° d = mo- The 
vectorial and axial interaction terms are neglected in Eq. (1). 

The quark fields are related to the gauge field A 11 through the covariant deriva- 
tive D^ = d* 1 - iA>*, where the gauge field is A^ = 5gA = -iS^A^ (the 
Polyakov calibration). The field $ is determined by tracing the Polyakov loop 

rP 



L(x) [5]: $L4] = ^-Tr c L(x), where L(x) = Pexp 



i I drA^x, t) 
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The gauge sector of the Lagrangian density (1) is described by an effective 
potential U (3>[A], $[A]; Tj fitted to lattice QCD simulation results in a pure 
SU(3) gauge theory at finite T [5,18] 
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Parameters of the effective potential (2) and (3) defined by fitting to the lattice 
results are summarized in [19] with the model parameter value of To = 0.19 
GeV. 

The grand potential for the PNJL theory in the mean-field approximation is 
given by the following equation [17] : 



fi($, <5, m, T,n)=U ($, $; T) + G(qq) 2 + Q q 



(4) 



with 
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E p is the quark energy, E p = y/p 2 + m 2 , E^ = E p =p fi (fi is the chemical 
potential) and N$(E P ) is the partition density with 



7V+(£ P ) = [l + 3 ($ + 5>e^) e"^ + + e" 3 ^ 
N$ (E p ) = [l + 3 ($ + $e _/3i ^) e"^ + e" 3/3 ^ 
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The gap equation for the constituent quark mass is obtained by solving the 
equation <9f2($, <§, m, T, /i)/dm = with the grand potential (4) 



m = m - N f G < qq >= m + 8GN c N f 



d 3 p 



m r 



/a (2tt) 3 E p 
where /$ , f$ are the modified Fermi functions 
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and /3 = 1/T is the inverse temperature. 
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Fig. 1. Meson mass M ff ff =t Y^ ff /2, the double quark mass and the meson width 
T n<a - in the PNJL model. 

The a and n meson masses are the solutions of the equation 



1 - 2G U ps/s (k 



0. 



where k 2 = m 2 and k 2 = m 2 in pseudo scalar and scalar sectors, and H ps / S 
are the standard mesonic correlation functions [20] 
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Both the pion-quark g^ qq (T, /i) and sigma-quark g aqq (T, p) coupling strengths 
can be obtained from U ps / S '- 



-2 rrp \ _ ^^■nqq/aqq\k ) , 

y-Kqqlaqqy 1 ' t 1 ) qj^ 2 |fe 2 =m| 
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The regularization parameter A, the quark current mass mo, the coupling 
strength G, the parameter To and the Mott temperature in such calculations 
are presented in Table 1. In the PNJL model there are two critical tempera- 
tures: The critical temperature of the chiral transition and the deconfmement 
temperature which can coincide for some model parameter set [19]. The first 
quantity is obtained as a minimum of dm/dT and the second one as a maxi- 
mum of <9$/<9T. For the vanishing chemical potential the temperature of the 
chiral transition is T c =0.213 GeV for our model parameters. 

The temperature dependence of n and a- meson masses is shown in Fig. 1, 
its behavior is typical for NJL-like models. The double mass of constituent 



Table 1 

The model parameters 



m [MeV] A [GeV] G [GeV]^ 2 T [GeV] T Mott [GeV] 



5.5 



0.639 



5.227 



0.19 



0.231 



quarks 2m is also plotted in this figure. As is seen, the curve for 2m crosses 
the 7r-meson one at T =0.231 GeV allowing the pion to dissociate into their 
constituents. This so-called Mott temperature is sometimes considered as a 
"soft" form of deconfinement. Note that Tm u > T c . At higher temperatures 
both scalar and pseudoscalar meson masses jointly increase. The decay width 
of ir and a mesons monotonically grows above the Mott temperature as T 
increases. 
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Fig. 2. Coupling constants g aqq (dashed line) and g ngg (solid line). 

As demonstrated in Fig. 2, in accordance with Eq. (11), the quark-cr coupling 
constant is rather weakly sensitive to the temperature slowly increasing with 
T but in the quark-pion case the coupling constant exhibits a kink singularity 
just at the Mott temperature. Technically, this feature results in the coupling 
strengths approaching zero for T — > Tyiott from below. This behavior differs 
markedly from the behavior of the couplings when evaluated in the chiral 
limit. 



3 Scattering cross section for qq — y qq and qq —y qq processes 



Quark-quark elastic scattering 



Let us consider now the quark scattering processes. The amplitude of the 
quark-quark scattering to 1/N C order is given by two diagrams presented in 
Fig. 3 taking into account channels with pion and cr-meson creation in an 
intermediate state: 
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where T n = {i^^-g^qq and r CT = l-g aqq . Here s and t are the usual Mandelstam 
variables. After appropriate transformations and bearing in mind that the 

total quark-quark amplitude is T qq = —— y^ \T t + T u \ 2 we get the following 

c c 

result : 



|T t | 2 =(|D?| 2 (*-4m 2 ) 2 + |D*| 2 t 2 ), 
|T B | a =(|I^| 3 (u-4m 2 ) 2 +raV), 

T<T " = "d^™*" + 4m2(M + t} " 16m2) " ^ CTjD > (t " 4m2) 
- DJD°t(u - 4m 2 ) + A^» , 

where the effective meson propagators in the t and u channels are 
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with coupling constants g( a ,ir)qq defined from Eq. (11) 
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Fig. 3. Diagrams of qq — > qq for ^-channel (left) and u-channel (right) 



The differential cross section has the form 

da P , \T\ 2 



dt 167rA(s,m 2 ) 



with A(s, m 2 ) = (s — 4m 2 )s. 

Then the total cross section of elastic scattering is: 



(T (A 



167rA(s,m 2 ) 



£ + dt \T\ 2 (1 - f F {£ T /i)) (1 - M^- T Ai)) , (19) 



where the inserted blocking factors fp take into account that rescattered par- 
ticles appear in a medium where other identical particles already exist. Here 
fp is the Fermi function fp = (1 + exp(/3x)) _1 and the sign in front of the 
chemical potential \x corresponds to a particle or an antiparticle. The inte- 
gration limits in Eq. (19) for the qq case with equal masses are t + = and 
t~ = Am 2 — s. The kinematic boundary reads s > 4m 2 . 

For calculation of the differential cross section in the center-of-mass system 
can be rewritten in the Mendelstam variables as t — —2p* 2 (l — cosO) and 
u = —2p* 2 (l + cosG), where in the center-of-mass system we have 



V = Qi = Q2 = 13 = ?4 



A 1 / 2 (s,m 2 ,m 2 ) 



(20) 



It is easy to show that the scattering amplitude for the process qq — >• qq can 
be obtained from the expressions for the process qq — > qq by the substitution 
t-H-t, s ■H- u , u ■H- s. The total elastic scattering amplitude for q q — > q q 
follows immediately from the qq — > qq amplitude by time reversal invariance, 
as long as no chemical potential is involved. 



Quark- antiquark elastic scattering 

Diagrams for the quark- antiquark process, qq — > qq, are given in Fig. 4. 
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Fig. 4. Diagrams of qq — > qq for t-channel (left) and s-channel (right) 



Similarly to the quark-quark case we have for the qq case 



|T t | 2 = (|Df | 2 (i - 4m 2 ) 2 + \Dl\H 2 ) , (21) 

|T S | 2 = (|D s CT | 2 (s - 4m 2 ) 2 + |Z>*|V) , (22) 



T t T s = -—(D^D^ts + 4m 2 (s + 1) - 16m 2 ) - D°D n s s(t - Am 2 ) 

- D?D° s t(s - 4m 2 ) + D*D* a ts) (23) 

with the meson propagators 

n 2 n 2 

1 ~ t-Mf ' s ~ s-M 2 ' [ ] 

(cr,7r) (cr,vr) 



Taking into consideration the isospin factors we can consider the following 
independent types of scattering reactions for quark-quark scattering : 



uu — > uu (dd — > dd), 

(25) 
ud — > ud (du — > du), 

where the first of them includes both channels and the second includes only the 
w-channel. In the same way, for quark-antiquark scattering we have a similar 
situation: 



uu — > uu (dd — > dd), 

ud — >■ ud (du^-du), (26) 

uu — > dd 

with both channels for the first reaction, the t- channel for the second one and 
the s-channel for the third one, respectively. 

The total and differential cross sections for antiquark-antiquark coincide with 
those for quark-quark scattering and can be calculated as in Eqs. (18), (19). 



4 Amplitudes for the qH — > qH process 



The Feynman diagrams for calculation of the quark-pion scattering amplitudes 
with the exchange of a quark, n- and u-meson are shown in Fig. 5. 

The amplitude in the s-channel (diagram 1 in Fig. 5) is as follows: 
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Fig. 5. The Feynman diagrams for s— and u— quark exchange channels (diagrams 1 
and 2) and the t-channel with the exchange by the sigma- meson (diagram 3). The 
double line corresponds to a pion. 
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After the transformation, Eq. (27) reads 



iT =- m (q 2 )^ 2 ( 275 ) (i7 B )^i«(qaJ 

s — m z 



(21 



where m is the quark mass. The result for the w-channel can be obtained by 
the substitution p± — > —p 2 
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The amplitude of the process with the a-meson exchange is as follows 



-iT 3 = u(q 2 )lu(q 2 )- —g a ^g aqq 



t- MV 



i T 3 = M fe) (70I70M92) 



1 



(30) 



t-Ml 



y<T7V7Ty(7qq t 



The total amplitude of the process is given as \T\ 2 = f c V] \T\ + T 2 + T 3 | 2 and 

c 

then after tracing and transforming equations we get: 
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T X T* = N c K 2 s [Mt - (a - m 2 ){u - m% (31) 

T x Tl=T 2 T{ = N C K S K U [-Mt + (s - m 2 )(u - m 2 )], (32) 

T t T; = T 3 T* = N c K s K t [m(s - u)}, (33) 

T 2 T* = N C K 2 U [M* - (s - m 2 )(u - m 2 )], (34) 

T 3 T; = T 2 T; = -N c K u K t [m{s - u)}, (35) 

T 3 T; = N c K 2 (4m 2 -t) , (36) 

where m, M n , M a are the masses of the quark, pion and sigma- meson, respec- 
tively. Summation over colors depends on reaction type. 

Here we have introduced the propagators K s , K u , K t 

2 2 

tv- _ 9-irqq ts- _ 9nqq ry- _ 9cnTiz9crqq / Q7 \ 

s ~ 2' u ~ 2' * ~~ 1 j\72~' y 6 ') 

s — m z u — m t — M£ 

where m, M n , M a are the masses of the quark, pion and sigma- meson, respec- 
tively. Here we have introduced the propagators K s , K u , K t 

2 2 

Ks = -^L_ # u = _?E££_ Kt = 9ann9 Z\ (38) 

s — m z u — m t — M z 

where g nqq , g aqq are defined by Eq. (11) and the coupling strength of a — nn is 
interrelated as g a7Tn = 2g aqq g 2 A aix -R and A a7Tn is the amplitude of the a — y tttt 
decay [21]. Summation over colors depend on the reaction type. Taking into 
account the diagrams with 1/N C and 1/N 2 corrections we have got color factors 
f c . Both color and flavor factors for every reaction are given in Table 2. 

Table 2 
Flavor and color structure of processes 



Process Isospin factor color factor f c 



U7T° — > U7T° 



\K S , \K u , K t [} + w c + lk 

U7T^ -> U7T~ K s , K u , K t I 1 + W 

^ d7r ° ^s,^ 2 K u ,K t (l + i 

d7r _^ dyr l Ks ^ l Ruj Kf (j + _2_ + ^ 

d ^+ __► dvr + K Ku , Kt ( 1+ i 



U7T 



iK s ±K u ,K t fl + i 



V2 s ' y/2 "' * V I N, 



Kinematic invariants for the qH — y qH scattering process are defined as fol- 
lows: 
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8= (qi + Pi) 2 = (g 2 +P2) 2 , 

* = (QI-Q2) 2 = (pi -P2) 2 , 

( V t v (39) 

u = {qi- P2) = (92 - Pi) , 
s + t + u = 2M 2 + 2m 2 . 



The differential cross section has the form similar to Eq. (18) 

d(Tel \T\ 2 



dt 16vrA(s,m 2 ,M2) 



(40) 



with A(s, m 2 , M 2 ) = (s — (M„. + m) 2 )(s — (M n — m) 2 ). Accordingly, the inte- 
grated cross section for the elastic qH scattering is: 

^ = 1fi w X 2 Tfis f dt l T | 2 ( X - f*( E * =F *0) (1 + /*(** T /*)) , (41) 

where i£ g , £"# correspond to the quark and hadron energies, the Bose-Einstein 
factor has the form /# = (exp(/3x) — 1) _1 and the integration limits are 

t ± = 2m 2 -^Us + m 2 - Ml) 2 T A(s, m 2 , M 2 )} . (42) 

As follows from Eq. (41) for the cross section, the g7r reaction has kinematic 
boundaries s > max {(m + M^) 2 , (m — M ff ) 2 }. 

For calculation of the differential cross section in the center-of-mass system 
we can use the same expressions as for the gg-scattering keeping in mind that 
p* 2 differs from Eq. (20) because the A(s,m 2 , M 2 ) factor has a different form. 



5 Numerical results and discussion 



Quark- (anti) quark elastic cross-section 



We shall start with the PNJL calculation of the in-medium qq, qq and qq cross 
sections. The integral cross sections o~ e i are plotted in Fig. 6 for all types of 
reactions shown in Eqs. (25), (26). 

The quark scattering cross sections are relatively featureless. For quark-quark 
scattering the flavor influences mainly the magnitude of the cross sections, 
being larger for quarks of the same flavors (a e i(uu — > uu) ~ 2.ha e i{ud — > ud)) 
(compare two upper panels in Fig. 6). This is because the ud scattering has in 
fact less exchange mesons available in the w-channel, since no neutral particles 
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Fig. 6. Total cross section of reactions qq — > qq (two upper plots) and qq —> qq (two 
bottom plots) at different temperatures. 

are admissible in this channel. The energy dependence of a e \ is very similar 
in both cases: the cross section has a clear maximum at the energy s/s ~1 
GeV for a moderate temperature T ~ 0.1 GeV which moves to smaller i/s as 
temperature increases exhibiting almost singular behavior at T = Tmou- The 
a e [ is getting rather flat if the temperature exceeds the Mott temperature (see 
the case T = 0.3 GeV in Fig. 6). 

Note that the cross section is evaluated according to Eq. (19,) where the Pauli 
factor for scattered quarks is involved. As is seen from Fig. 6, a e \ decreases 
when y/s grows which is a consequence of nonperturbatively in the coupling 
constant G of the PNJL treatment. In the Born approximation the inverse 
behavior a e i ~ s is observed, as was demonstrated in Ref. [23]. The behavior 
of the scattering cross section for quark- ant iquark of different types is very 
close to that for quark-quark of different flavors (compare ud and ud reac- 
tions in Fig. 6). However, for quark- ant iquark of the same flavor (uu and dd) 
the cross section at any temperature has the resonance-dominated behavior 
demonstrating a huge maximum located at the energy close to the u-meson 
mass. In other words, the quark- ant iquark scattering shows a threshold di- 
vergence at the Mott temperature Tmou-, a t which the pion dissociates into 
its constituents and becomes a resonant state. This feature manifests itself in 
other processes like qq — > 77 [25] as well as tttt — > nir [21,26] and tt'j — > 717 [27]. 
Dramatically high cross sections mean that near Tmou & local equilibrium can 
be established. 



13 



s 



_.w 














uu->uu 










1.5 


SU/(3) _..- 








..--' '. 












1.0 


" 








■ 




-•'' 




PNJL, 


T: 


= 0.27 GeV : 


0.5 




n n 





SU/(2) 




. i . 








0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Vs" (GeV) 

Fig. 7. Comparison of the energy dependence of the elastic cross section for the 
uu — > uu scattering for a temperature above the Mott temperature. The dotted 
and dashed lines are the NJL results for SUf(3) and SUf(2), respectively, from 
Refs. [24,23], the solid line is our two-flavor PNJL calculations for T = 0.27 GeV. 

The calculated results can be compared with those obtained earlier [24,23]. 
Some elastic scattering cross sections were calculated in the two-flavor sector 
SUf(2) of the NJL model with the additional restriction of the chiral limit 
condition m =0 in Ref. [23]. The three flavor NJL model was considered 
in [24] and the results are exemplified in Fig. 7. R is seen that the SUf(3) 
calculation yields a larger cross section for uu ~ uu than the corresponding 
two-flavor case. This is due to the point that the additional exchange channels 
t], rj' and a' are missing in the two flavor case. In particular, the uu — > uu 
channel of the elastic scattering at T = 0.215 GeV for SUf(3) is regularly 
above that for SUf(2) by the factor of (3-4) in the whole energy range from 
the threshold till yfs ~1.2 GeV [24] (cp. two dotted lines in Fig. 7). The values 
of the uu — > uu cross section for SUf(2) [24] with Tmou are consistent with our 
PNJL results presented in Fig. 7 if one takes into account the different Mott 
temperatures in these calculations. R is not surprising since the suppression 
due to the Polyakov loop works only in the low temperature sector, T < Tu tt- 

For each process the differential cross sections may also be calculated. As 
an example, the angular distribution in the elastic qq and qq scattering is 
presented in Fig. 8 at the energy -y/i =1 GeV. In the case of quarks of the same 
flavors (uu), the angular distribution is isotropic except for the narrow region 
|cos0| ^0.9. The differential distributions for quarks of various flavors (ud) or 
different types (ud) are also almost flat besides the region |cos0| >0.8 but this 
deflection from isotropy is in backward direction for ud while it is in forward 
direction for ud reactions. For uu and dd scattering, the angular distribution 
changes from isotropic to clearly anisotropic one if the temperature decreases 
from Tmou to about 0.1 GeV. In accordance with the temperature dependence 
of a e [, the magnitude of the differential distribution of da e i/dcosQ decreases 
with the growth of the temperature. 
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Fig. 8. Angular distribution of the elastic scattering processes qq — > qq (upper 
panels) and qq — > qq (bottom panels) at ^fs = 1 GeV for different temperatures. 

The processes presented in Figs. 3 and 4 can also be calculated in the lowest 
order perturbative QCD. These on-shell high-energy calculations are in use 
starting from the first classical parton cascade model [28]. In Fig. 9, the com- 
parison between the PNJL and the pQCD results (obtained according to the 
model description in Appendix of Ref. [23]) is given. As is seen, in contrast 
with the PNJL, the pQCD description predicts a strong scattering enhance- 
ment at forward angles for finite parton masses. In the limit m — >0 the cross 
section increases and diverges. 
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Fig. 9. Comparison of the differential ud —> ud cross section calculated in the pQCD 
model at y/s =1 GeV for finite (dashed line) and vanishing (dotted line) masses. 
Our results for T = 0.25 GeV are plotted by the solid line. 

Such behavior is more transparent from simplified two-body pQCD calcula- 
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tions with m = used in the AMPT model [9] . In particular, for gluon elastic 
scattering (which differs by the Casimir factor from qq scattering), one has 

d<jgg Qua 2 ut us st Qna 2 1 1 Qna 2 s 

with the strong coupling constant a s . Equation (43) is obtained in the leading- 
order QCD by keeping only the leading divergent terms for identical particles, 
which allows one to limit oneself to the angle range < O < it/2 (the last 
equality in (43)). The cross section really diverges at the scattering angle 
O =0. The singularity in this cross section can be regularized by introducing 
the Debye screening mass Mp, leading to 



d(J 99 = 9to s (u] 

dt 2(t-M D ) 2 K ] 

and respectively for the total cross section at relativistic energies 



9ira 2 1 9-na 2 

a " ~ 2M 2 D 1 + Mljs ^ 2M| ' 

if s ^> M D [9]. Note that in this approximation the pQCD gives energy- 
independent cross section (45) in disagreement with the NJL-like chiral model 
discussed. Taking Mr, =3 fm _1 we get a e \ =3 mb. In the real AMPT calcu- 
lations Mo is a parameter and this cross section changes from 3 to 10 mb in 
different model versions [91. 



Quark-hadron elastic cross-section 



The integral cross section for the in-medium elastic scattering of quarks on 
pions is presented in Fig. 10. For the considered reactions un~ , dir + , un and 
their charge conjugated the cross sections behave very similarly: a e i is maximal 
just at the threshold Sthr — (m + m^lT)) 2 reaching here values of 60-80 mb 
and then monotonically falls down rather quickly with the slightly different 
slopes for different reactions. At the Mott temperature and for all energies the 
qH elastic cross section exhibits a huge maximum as large as several hundred 
of mb (to be cut in Fig. 10). Generally, the magnitude of the quark-pion cross 
sections is higher than the quark-quark one and comparable with free hadron- 
hadron cross sections. In the case of the uir~ reaction, the cross section fall-off 
is the slowest and the quark-pion scattering turns out to be by several times 
higher than that for quark-quark in the energy range ~(0.7-1.5) GeV. 

As is seen from Fig. 10 (right-bottom panel), o~ e i(T) is weakly changing func- 
tion up to T ~ 0.18 GeV and then has a kick at the Mott temperature. 
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Fig. 10. Energy dependence of the total scattering cross sections for the processes 
qH — > qH, corresponding to the different reactions and the temperature dependence 
of different components (see the text) of the <i7r° — > uir~ scattering cross section at 
y/s =1 GeV (right bottom panel). 



This cross section is the sum of different terms corresponding to the squared 
amplitudes: an (see Eq. (31)), 022 (Eq. (34)), 033 (Eq. (36)) and summary in- 
terference term \ nter / (Eqs. (32) + (33) + (35)). The dominant u-channel gives 
CT22 ~ 0.5 a e \. In contrast, the cross section of the t-channel is very small, 
0-33 —0.5 mb. The rest of a e \ is shared almost equally between the sum of 
all interference terms (J inter f (note the scale in figure) and the s— channel o"n 
elastic cross section. 

In Fig. 11, some examples of the angular distribution are given for the quark- 
pion scattering. At the energy yfs < 0.5 GeV independently of T besides the 
vicinity of Tmou, the angular distributions is isotropic with high accuracy. At 
higher energy and moderate temperatures T < Tm u the differential distri- 
butions are almost isotropic with small anisotropy growing at the backward 
scattering angles. This backward bounce is more prominent for lower energies. 
The presented angular distributions really are the contributions of diagrams 
of three types (see Fig. 5). As is shown in the right-bottom panel of Fig. 11 
and in accordance with the above discussion, <r 2 2 is dominant and completely 
defines the backward scattering (note the scale factor 0.2 in the figure). The 
aii(G) and 033(0) angular distributions are forward peaked but the asymme- 
try is much more pronounced in the s-channel. The sum of interference terms 
results in the distribution with a broad maximum at cosQ ~0. 
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Fig. 11. Angular distributions of the uir~ — > wr~ scattering for all reactions at 
various temperatures and energies and the contributions of different channels of the 
d,7r° — > wk~ reaction to da e i/dcos@ for sfs =1 GeV and T = 0.1 GeV (right-bottom 
panel). Notation of components is the same as in Fig. 10. 

The presented results are the same for the uix~ — > wk~ and 7r + — >■ dir + reac- 
tions because their flavor factors are the same (see Table II). For other channels 
of the quark-pion scattering, the shape of angular distributions is similar to 
that for U7i~ — > wrr~ but the magnitude is lower in accordance with the flavor 
factors in the Table II. Because a pion consists of a quark-antiquark pair, the 
substitution u — > d for quarks does not change the calculated scattering cross 
sections results in Fig. 10 and angular distributions in the s- and u-channels 
shown in Fig. 11 are exchanged. 



6 Concluding remarks 



In- medium elastic scattering of quarks is evaluated by treating quarks in terms 
of the two-flavor PNJL model accounting for their scalar and pseudoscalar 
interactions. The model parameters are in agreement with low-energy static 
mesonic properties and lattice QCD results for gluons. General trends of the 
energy and temperature dependence of a e i and da e i/dcosQ are investigated 
in a large range of s/s and T at both below and above the Mott temperature. 
Essential influence of the flavor of interacting quarks on the scattering results 
is noted. The comparison with the earlier works shows that the SUf(2) PNJL 
results agree with that of the NJL model for the same SUf{2) symmetry (at 



T > Tmou) but noticeably below in the case of three flavor. This difference 
is due to a larger number of contributing intermediate states for the quark 
scattering in the case of SUf(3) symmetry. 

The above PNJL results are presented for the vanishing chemical potential. 
However, \i enters not only into the Pauli exclusion factor but also into the 
polarization function Il(/c 2 ) which makes masses and coupling constants to be 
/i-dependent. To consider the quark-baryon case, an additional issue arises: one 
should treat properly the loop formed by the quark and the diquark that form 
the baryon. As demonstrated in [29], taking into account the finite chemical 
potential results in some suppression of the scattering cross sections at large \i. 
The PNJL approach should also be generalized to the SUf(3) symmetry. This 
allows one to extend the set of reactions including strange quark and strange 
hadrons and to effect non-strange reactions due to increasing a number of 
possible intermediate states, as noted above. This work is in progress now. 

First predictions for the quark-pion scattering are given. The integral cross 
sections for this reaction are higher than those in the quark-quark case. Quark 
differential cos0-distribution is practically constant with some enhancement in 
the backward direction at ^Js >1 GeV. The developed technique and obtained 
results can be applied in the kinetic approach [9,11,14] to take into account the 
interaction between quark-gluon and hadronic phases. The derived formulae 
for elastic quark-hadron scattering can be easily generalized to the expression 
for gamma and dielectron emission [30]. This new channel can definitely give 
rise to an observable effect and would confirm the existence of the quark- 
hadron interaction. 

Another possible implementation of these results is a study of transport prop- 
erties of the system. For example, the shear viscosity 7] can be estimated in the 
so-called relaxation time approximation by using the average momentum loss 
p, the quark densities rij and the mean life time as rj ~ V] riipi/Xi, the mean life 

i 

time being inversely proportional to the quark cross sections A -1 = y^ nja Q j 

i 

Ref. [31]. In the two-flavor NJL model, this shear viscosity was estimated in 
Ref. [23] for the chiral limit m = 0. One should not expect a noticeable dif- 
ference with our model treating a quark system with the finite quark masses 
in the SUf(2) PNJL model, especially at high temperature where the NJL 
model is more justified. However, for the mixed quark-hadron system consid- 
ered in [9,11,14] the shear viscosity should be smaller, because the qH elastic 
cross section is noticeably larger than quark-quark one and particle density of 
quarks and pions is rather abounded in the mixed phase. Certainly, this prob- 
lem deserves a more elaborated study, for example, by the method developed 
in Ref. [32] . 
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